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Abstract 

We present a study of hopping conductivity for a system of sites which can 
be occupied by more than one electron. At a moderate on-site Coulomb 
repulsion, the coexistence of sites with occupation numbers 0, 1, and 2 results 
in an exponential dependence of the Mott conductivity upon Zeeman splitting 
fisH- We show that the conductivity behaves as In a = {T/To)^/^F{x), where 
-F is a universal scaling function of x = fiBH/T{To/T)^/^ . We find F{x) 
analytically at weak fields, x <C 1, using a perturbative approach. Above some 
threshold Xth, the function F{x) attains a constant value, which is also found 
analytically. The full shape of the scaling function is determined numerically, 
from a simulation of the corresponding "two color" dimensionless percolation 
problem. In addition, we develop an approximate method which enables us 
to solve this percolation problem analytically at any magnetic field. This 
method gives a satisfactory extrapolation of the function F{x) between its 
two limiting forms. 

PACS Numbers: 73.40, 71.25 
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I. INTRODUCTION 



Low-temperature conductivity in a disordered semiconductor is controlled by phonon- 
assisted electron hops between localized states. At sufficiently low temperatures, only those 
sites which have energy levels close to the Fermi level participate in the hopping transport. 
This defines the Mott variable range hopping ( VRH) regime [|[] . Mott conductivity depends 
exponentially on temperature, (j{T) oc exp [-(To/T)^/('^+i)], where Tq = pd/ga'^ is the char- 
acteristic Mott temperature, g is the density of localized states at the Fermi level, a is the 
localization radius for a single site, d is the dimensionality of the sample, and the numeri- 
cal factor jSd is determined by percolation theory 0. In the more standard case of lightly 
doped semiconductors, the strip of localized states in energy space is relatively narrow, and 
each site can accommodate at most one electron. Under these conditions, the spin degree 
of freedom of the hopping electrons has no effect upon the exponential factor of the hop- 
ping conductivity. The application of a magnetic field H results in the modification of this 
exponential factor due solely to the orbital effect of the field 0. 

Kamimura et al. were the first to recognize that the spin degree of freedom plays a 
significant role in the magnetoresistance if a certain fraction of the sites can accommodate 
more than one electron. Double occupancy is possible if the on-site Coulomb repulsion U 
between the electrons is smaller than the width of the distribution function of the energies 
of the localized sites. In this case there are two types of sites which contribute to hopping 
transport. Sites of the first type, which we will call type A, have energies e close to /i. Sites 
of type B have one electron at a deep level with e ~ — ?7, so that the energy for the second 
electron is close to /z. The sites which are neither of type A nor of type B have energy levels 
which are too far from the Fermi level to contribute to transport. At zero magnetic field, the 
probability for two electrons on two singly-occupied sites to have opposite spins equals 1/2. 
Therefore, hops between A and B sites can occur. In the strong field limit however, all spins 
are polarized, and A ^ B hops are completely suppressed assuming that two electrons 
occupying the same site form a singlet state at all relevant magnetic fields [Q. Thus, at 
sufficiently strong magnetic fields the characteristic Mott temperature is determined by the 
larger of the two densities of states gA and gB rather than by the net density of states 
g = gA + gs- The increase in Tq due to the field induced suppression of A B hops leads 
to a giant positive magnetoresistance. 

In this paper we present a detailed study of the hopping magnetoresistance induced by 
Zeeman splitting. We show that the criterion for the strong magnetic field limit described 
above corresponds to a finite threshold value Ht^. Below this threshold value, the conven- 
tional Mott exponent, {Tq/TY^^^'^^\ is modified by a factor F{x), 



The universality of the Zeeman splitting induced magnetoresistance is a key result of this 
paper. Its significance is illuminated by noting that it reduces the calculation of the mag- 
netoresistance to the determination of a universal function F{x) of a single dimensionless 
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which is a universal function of a single scaling parameter 
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parameter x, as opposed to previous solutions 0] for which the magnetoresistance was cal- 
culated function of two variables, T and H. 

Above a certain threshold x > Xtt, the universal function F{x) attains a constant value 
which we determine analytically. In addition to the saturation value, we find the x -C 1 
asymptote of F{x) analytically using a perturbative approach. In an attempt to obtain an 
analytically determined fitting curve for the universal scaling function F{x), we extend the 
invariant technique to the case of the "two-color" percolation problem. To determine the 
accuracy of these two analytical techniques, we find the scaling function F{x) numerically by 
simulations of the corresponding "two color" dimensionless percolation problem and compare 
it with the functions obtained analytically. 

Some of the results of this study have been reported in a short communication P|. 
Another mechanism of giant magnetoresistance due to two types of localized states was 
considered in Ref. |0. 



II. APPLICABILITY OF THE MODEL 

Zeeman splitting can make the dominant contribution to the magnetoresistance in a 
number of important cases, which include those of undoped amorphous silicon and 
moderately disordered two-dimensional electron systems. In the latter case, the minima 
of the random potential may serve as sites accommodating several electrons. The orbital 
effect can be avoided altogether if the field H is applied parallel to the plane of the two- 
dimensional system. In amorphous silicon electrons are localized at "dangling bonds" [§. 
For the case of bulk a-Si, the orbital effects of the magnetic field cannot be eliminated, 
and Zeeman splitting makes the dominant contribution to the magnetoresistance only in a 
certain range of temperatures. To determine this range, we consider the tunneling of an 
electron a distance L through a barrier consisting of the superposition of an "intrinsic" , Vq, 
and a "magnetic", ma;^a;^/2, barrier {ujc = eH/mc is the cyclotron frequency of the field 
H). The action for the electron subbarrier motion is 



- = -/^2™(t/„ + -™.J.^)d.«- + — . (3) 
The localization radius of the electron in the absence of a magnetic field is related to the 



strength of the "intrinsic" barrier Vq as a = h/ \/2mVo, and A = ^JcTi/eH is the magnetic 
length. The second term in Eq. (0) represents the correction to the subbarrier action due to 
the presence of a magnetic field. From Eq. (^), we see that the orbital effect of the magnetic 
field is neghgible at 

The characteristic tunneling distance L in the VRH regime is of course determined by the 
temperature 

, i/(d+i) 

L ^ a (^) . (5) 
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Eqs. (D and (^) determine the upper limit of the magnetic field. Requiring the Zeeman 
splitting to be strong limits the magnetic field from below: 

T « ^bH « ^ O-) . (6) 



2ma'^ \ToJ 
This relation places the restriction 

_2 X 2(d+l)/(2d-l) 



T<To[ 1 . (7) 



upon the range of temperatures at which Zeeman splitting makes the dominant contribution 
to the magnetoresistance. 



III. DERIVATION OF THE CONNECTIVITY CONDITIONS 

Mott's arguments enable one to find the temperature dependence of the exponential 
factor in VRH conductivity. However, it is necessary to use percolation theory to determine 
the numerical factor Pd in the exponent Similarly, to study quantitatively the effect 
of Zeeman splitting upon the conductivity one has to reformulate the problem in terms 
of percolation theory. Because two different site types are now involved, the percolation 
network which determines the hopping conductivity consists of three types of links: AA, BB, 
and AB. In order to find the VRH conductivity, we must first determine the conductances 
Gaa, Gbb, and Gab of the three types of elementary links. To begin with, we introduce the 
probabilities Pi of having the site i occupied by zero (0), one (t or |), or two (t|) electrons: 

Pi(0) = Z-\ (8) 

P,(T) = Z-'exp[-{e,+ftBH)/T], (9) 

P,(i) = exp[-(e, - iibH)/T], (10) 

mi) = ^-iexp[-(2e,, + U)/T]. (11) 

Here 

Z = l + exp[-(e, + fisH)] + exp[-(e, - fisH)] + exp[-(2ei + U)/T]. (12) 

Neglecting the pre-exponential factors (which include, e.g., the deformation potential con- 
stant), we can express the three elementary link conductances as products of phonon and 
electron occupation numbers. 

Gaa oc [Pi(T) + Pi(i)]P2(0)iV(e2 - e,) exp(-2P/a), (13) 
Gbb oc Pi(n)[P2(T) + P2(i)]iV(e2 - d) exp(-2P/a), (14) 
Gab oc [Pi(T)P2(i) + Pi(i)P2(T)]iV(e2 + U ~ e,) exp(-2P/a). (15) 

Here R is the distance between two sites and is the Bose distribution function. 

The standard approach to the formulation of the percolation problem requires the ex- 
ponential representation of the elementary link conductance, G oc e~^. In this colored 
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percolation problem, there are three different exponents ^aa, ^bb, and ^ab which are ex- 
tracted from Eqs. (|13|)-(|15|). Assuming that the on-site Coulomb interaction is very strong 
U > iibH > T, we find 



^AA 
C.BB 



1^2 - eil + |6l - flBH\ + \€2 - ij,bH\ W. 

2T ^ a ' 

|e2 - eil + |ei + f/ + fiBH\ + \e2 + U + ^ibH\ 



2T 



+ - 



2R 



lAB 



\e2 + U -ei\ + |ei - iibH] + \e2 + U + ^ibH\ 



2T 



HbH 2R 
T a 



(16) 



(17) 



(18) 



The energies ei and €2 in Eq. (|TB[) correspond to the energy levels of sites of type A and B 
respectively. The Fermi level depends upon magnetic field due to transitions from doubly 
to singly occupied sites, 5/i = ^bH^Qb — gA)/{gA + Ob)- However, this has no effect upon 
the resistance, provided that qa and qb are energy independent for energies on the order 
of fiBH- Unlike the case of semiconductors doped by shallow impurities [|], in amorphous 
silicon the densities of states are independent of energy on this scale, since they vary with 
energies comparable to the onsite Coulomb interaction U ~ 100 meV 0. Therefore, we can 
safely make the assumption of constant densities of states in the vicinity of the Fermi level. 

Instead of the energies it is convenient to introduce a new set of variables Ei, chosen 
in such a way that their values are close to zero for sites participating in electron transport: 



We further define 



e = e — iibH for A-sites, 

e = e + U + i^bH for 5-sites. 



,TT\ \e2- ei-2^.BH\ + \ei\ + \e2\ , „ 
ei2[H) = + hbH- 



(19) 



(20) 



The new variables simplify the dimensionless exponents (|T( 

£12(0) 



^AA — t 



BB 



T 



2R 
a 



Wi) to the forms: 



(21) 



6i2{H) 2R 

Ub = ^^ + —. 22 
I a 

In the absence of a magnetic field, C,aa = (,bb = (,ab, and we return to the standard 
percolation problem for VRH conductivity. At if = only the net density of states g = 
gA + Qb at the Fermi level is relevant, and the ratio of the densities of states Qa/qb does not 
affect the conductivity. 
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IV. MAGNETORESISTANCE IN THE LIMITS OF WEAK AND STRONG 

FIELDS 



A. Strong Field Limit 



As follows from Eq. (^), the exponent ^ab can not be smaller than 2^bH/T . Thus at 
sufficiently strong magnetic fields the transitions between sites A and B cannot occur In 
this limit, the conductivity is determined by two parallel percolation networks, one of which 
consists only of type A sites and the other of only type B sites. Therefore the conductivity 
is independent of H and satisfies Mott's law with the density of states g = ma.x{gA, Qb}- 
This picture is valid above a certain threshold for the magnetic field: 



T ~ 
H > Hth = - — 



(23) 



when all possible C,ab exceed the critical exponent determined by the percolation theory 
solution of the VRH problem [Q], 



{gA + 9b)To 
~gT 



i/id+i) 



/3d 

(qa + gB)a'' 



(24) 



B. Weak Field Limit 



The opposite limit of weak fields also allows for analytical consideration by means of 
a perturbative approach applied to the percolation problem with the density of states g = 
QA + gB and the percolation threshold Q = (To/T)^/('^+^\ At if -C Hth possible field induced 
variations of ^ab are small, A^^^ -C This enables us to use the perturbative approach 
proposed in 0. According to Ref. 0, one can find the small shift of the percolation threshold 
as an average increment of caused by the small perturbation, A,^c = {^0 ■ To calculate 
the average (A^) over the statistical ensemble of sites, one needs to find the correction A^ to 
the exponent for each link due to a small magnetic field. Clearly, A^aa = A^bb = 0, whereas 
A^AB = '^fJ'BH/T if El > 62, and A^^^ = otherwise, see Eqs. ( pi| ) and (^). Taking into 
account the fraction of AB links with ei > 62 in the percolation cluster gAgB/{gA + gB^, 
we find 

The dependence of the hopping conductivity upon the critical exponent is cr ~ e'^". 
Therefore Eq. (p5[) enables us to find the dependence of the conductivity on the magnetic 
field: 

a{T,H) 2gAgB I^bH . 

"'7m=-(5I^— '''' 



This dependence is sensitive to the relative densities of A- and -B-sites. As expected, there 
is no field dependence at = or = 0. 
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The region of validity of ( p6D is determined by the apphcabihty of the perturbative 
approach, i.e., by the requirement A^^^ <^ In terms of the magnetic field strength, 
this condition reads fisH <^ /is-f^th ~ TC,^. Since ^ 1, the latter condition does not 
contradict our initial assumption, that fisH 3> T. More precise limits of applicability of the 
perturbative calculation become apparent from the comparison of Eq. p6| ) with the results 
of numerical simulation which we present later in this paper. 



V. SCALING ANALYSIS 
A. Scaling Conjecture 

It is worth noting that the dimensionless parameter of the perturbation theory is 

The threshold field i^th corresponds to the universal (i.e., temperature independent) value 
of this parameter, 

= 2 [-^) ■ 

The existence of the dimensionless parameter x allowed us to make the conjecture (^, 
which corresponds to the following scaling behavior of the percolation threshold: ^c{T, H) = 
^^F{x). This scaling function F{x) determines the conductivity at finite magnetic fields 

From the cases of low and high fields discussed above we already know the limiting behavior 
of the function F{x). It has linear expansion at small x and reaches a constant value at 
large x: 

f 1 + T^^r^a; at x < 1, 



B. Proof of the Scaling Conjecture 

We will now prove our conjecture by reducing the initial problem of hopping 

magnetoresistance to a dimensionless percolation problem. 

We start from the conventional percolation approach to the hopping conductivity and 
introduce a positive variable ^. At given ^, all links are cut except for those with conductance 
G > e~^. At small ^ percolation does not occur, but at some particular value ^ = the 
network starts to percolate. This threshold value determines the conductivity a ~ e'^". Our 
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goal is to find the dependence of on T and H. We generalize the approach of Ref. and 
introduce a set of dimensionless variables A, p, and x defined by the following relations: 

e = T^A, R=^a^p, pBH = Tix- (30) 

In these variables the connectivity condition for a link is: 

Ai2(0) + p < 1 ioi AAk BB links, . . 

Ai2(x) + P < 1 ioi AB links ^ ' 



where 



|Ai| + |A2| + |Ai-A2-2x| , 



Clearly, only sites with the dimensionless energies |A| < 1 can be connected. The dimen- 
sionless concentration of these sites in p-space is 

n = ^^{gA + gB)a''Te+\ (32) 

We are now in a position to formulate the dimensionless percolation problem associated 
with the hopping magnetoresistance problem. Consider a random distribution of points 
with concentration n in a dimensional p-space. Each point is characterized by its type 
[A or 5), radius- vector p^, and energy Aj. The latter is distributed uniformly over the 
interval (—1,1). Two points form a link if the condition (|3lD is satisfied. The problem 
is characterized by two dimensionless parameters: the ratio 7 of concentrations of A- and 
5-points, and dimensionless magnetic field x- At 7 = Qa/ Qb this dimensionless problem is 
equivalent to the original percolation problem for the hopping magnetoresistance. To reach 
the percolation threshold, we increase the total dimensionless concentration n (holding 7 
constant) until the critical value nc(x) is reached. Once ^^(x) is found, we can determine 
the threshold for the original problem from the relation (p2D, which can be rewritten as 

nc(x)K(0) = {UilY^'. (33) 

Rewriting x as defined by Eq. (|30| ) in terms of the parameter x (see Eq. (|27|)) we find that 
must be a solution of the following equation: 



(34) 



One can easily show that this equation has exactly one solution. Then it obviously has the 
form 

^c = ecF{x). (35) 
This proves our scaling conjecture (|I]), (||). 
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VI. SIMULATION 



As was previously mentioned, the accuracy of the perturbative approach used at small 
X, see Eq. ( ]29| ) can only be determined by comparison with simulation data. In addition, 
there is an intermediate range of field strengths within which F{x) must be numerically 
determined, as the analytical theory does not extend to this range. To perform these tasks, 
we developed a code which closely resembles that of Skal et al. [0. We will describe the 
method as applied in two dimensions. The simulation consists of the following: N sites are 
randomly thrown in an / x / box (which is large enough for the system to be well below the 
percolation threshold for the given value of A^). Each of these sites is randomly assigned a 
site type, A oi and an energy A. Two sites are connected if their parameters satisfy the 
condition (0). Percolation is said to occur when two strips (one at each end of the box), of 
width equal to the average distance between sites, are connected. We start with / 3> lc{x)i 
and make I smaller until this connection occurs. This determines the percolation threshold 
Icix)- The critical site concentration is ndx) = ^/^lix)- 

We ran the simulation at 7 = 1 and at 7 = 1/2 on 3600 sites and averaged each point over 
100 runs. One can see from the data shown in Fig. 1, that F{x) is an increasing monotonic 
function which attains a limiting value for x > Xth- At 7 = 1 the limiting value of F{x) 
is systematically suppressed below the exact result given by Eq. (P^). This is an artifact 
of the simulation which occurs because 7=1 corresponds to the same number of A and B 
sites. For x > Xth at 7 = 1, there are two parallel networks which can percolate. Thus for 
X > Xth, the simulation is in effect being run twice, once on the A-network and once on the 
i?-network, with the smaller value of the two resulting concentrations selected. Clearly, this 
finite size effect does not exist at any other value of 7. The simulation data of Fig. 1 used 
in conjunction with the scaling functional form (|l|) yields the magnetoresistance for 7 = 1 
and 7 = 1/2 at arbitrary fields. 

In order to determine the range of our small x asymptotic form, we replotted the data 
as critical site concentration ric vs. x ^-nd found that nc{x) begins to deviate from being 
linear in x fo^' X ~ 0.20. Weighted linear fits to the data yielded the zero field values 
n^(0) = 7.064 ± 0.023 at 7 = 1 and 71^0) = 7.036 ± 0.031 at 7 = 1/2 both of which are 
in good agreement with the accepted value of 6.9 ± 0.4 |1^. To test the accuracy of our 
weak field perturbative approach, we performed a weighted linear fit to the data taken at 
and below x = 0.20. As one can see from the insets in Fig. 1(a) and 1(b), the fits to the 
simulation data yielded slopes and intercepts which were in excellent agreement with those 
obtained perturbatively at both 7 = 1 and 7 = 1/2. 

VII. APPROXIMATE INVARIANT FOR THE PERCOLATION PROBLEM 

Ideally we would like to know the scaling function F{x) for any value of 7. However, the 
simulations used to determine F{x) for each value of 7 require a large amount of computer 
time. Thus we only found F{x) for two values of gamma: 7 = 1 and 7 = 1/2. The devel- 
opment of an analytical method which yields an approximate form of the scaling function 
F{x) is therefore highly desirable. The invariant method developed by Shklovskii and Efros 
gives just such a simple approximate solution of the VRH problem. Unfortunately, this 
approach cannot be adapted to the case of the two color percolation problem. In this paper 
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we suggest an alternative invariant. This approach (i) gives a better estimate of the constant 
for the standard percolation problem, and (ii) can be generalized to the two-color case. 
As it reduces the problem of determining the universal scaling function F{x) to the solution 
of a simple integral equation, we have used it to obtain an analytically determined fitting 
curve for the universal scaling function F{x). 



Recently the same invariant conjecture was independently proposed by loselevich 11 



A. Invariant approach to the standard VRH Problem 

We will first formulate our invariant conjecture as applied to the standard (single-color, 
93 = 0) VRH problem. At the percolation threshold (^ = ^c), the probability density of 
finding a defect state with energy 62 linked to one with energy ei is 

GAAiei,e2) = J 0{ic-UA{ei,e2,r))d'^r 

.,„,ff|Vfi_£HMyv«.--MV (36) 



Here Vd = 27^'^^'^ /T(d/ 2) d is the volume of a d-dimensional sphere of unit radius. For the VRH 
problem, the energies of the two sites are not fixed but rather are governed by distribution 
functions. Therefore the normalized distribution of energies fi{£2) of sites which can be 
connected to a site in the system is 



^1/1(62) = J &iGaa(£i,£2)/o(^^i) 




9AV, (^) / (1-^M(1- ^) fo{e^)de,. (37) 

Here /o(ei) is the distribution of energies of the first site; the normalization coefficient 
Ti is the average number of sites which can be directly linked to the first one (i.e., the 
average number of bonds per site). Repeating the procedure of Eq. (|37|) several times, we 
find that the distributions of site energies obtained in this manner converge rapidly to the 
eigenfunction of the dimensionless integral equation 

A,/(A2) = /[I - Ai2(0)]"^^(l - Ai2(0))/(Ai)rfAi. (38) 

We conjecture that this eigenfunction is approximately equal to the distribution of energies 
for the sites in the infinite cluster. The eigenvalue is proportional to T which would then 
be the average number of bonds per site in the infinite cluster, 

T = Xd9AVd(^^^'T^c. (39) 

We further conjecture that T is an invariant of the percolation problem, and equals the 
average number of bonds per site Be for the random sites (RS) problem [0. With this 
assumption, we obtain 
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where and nj are the critical concentrations of the RS and VRH percolation problems. 

A numerical solution of the integral equation (|38| ) gives: A2 = 0.4301, A3 = 0.3154, A4 = 
0.2489, in excellent agreement with the values given by Eq. (0), with nf^ and n^^^ obtained 
via simulations of the corresponding percolation problems: A2 = 0.410 ±0.004, A3 = 0.303 ± 
0.004, A4 = 0.232 ± 0.005. The simulation procedures we used to find n^^ and n^^^ are 
nearly identical to the one described in Section and were run on = 8100 sites and 
averaged over 50 runs for d = 3 and d = 4 and on = 6400 and averaged over 100 runs for 
d = 2. 



B. Generalization to the "Colored" Percolation Problem 

We shall now generalize the discussion of the previous section to the case of a two-color 
model, qb ^ 0- The sites are now characterized by two parameters: a continuous variable 
e and a discrete variable, A or B. Consequently, instead of a single equation (|38|), we now 
have a system of two integral equations: 

T ( {^\'^\ ]= [( G^^f^f) GAB{e,,e,) \ / Us.) \ , (4,) 
V/b(^2) J J \1 GAB{e2,ei) 7 GAA{ei,e2) J V / 

Here, Gaa{£i, £2) is given by Eq. (P^). The off-diagonal element Gab{£i, £2) is the probability 
density of finding a link between an A-state with energy ei and a 5-state with energy £2- 
It can be obtained from the probability density Gaa given by Eq. ( ^61) by making the 
replacement £12(0) EuiH). 

The kernel of this integral equation is not symmetric. However it can be symmetrized 
by performing a linear transformation: 

(42) 

Clearly, this transformation preserves the eigenvalues. In the new basis, the integral equation 
(p]) attains a symmetric form 






G'aa(£^1,£^2) 7 ^^'^GABi£l,£2 

7"^^^G'ab(£2,£i) l^^GAAiSi,e2) 




(43) 



guaranteeing that all of its eigenvalues are real. 

We now extend the method of invariance to the "two-color" percolation problem and 
conjecture that the maximal eigenvalue T is independent of the magnetic field. Let us show 
that by assuming the invariance of T we recover the results of the weak field perturbative 
approach described in Section |1V Ij| . First we solve equation ( ^3D at H = 0, for which 
Gab{£i,£2) = Gaa{£i,£2)- One can easily see that the normalized solution corresponding 
to the maximal eigenvalue has the form 
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Here the function f{e) is defined as the normahzed solution of the "monochromatic" integral 
equation 

Ao/(e2) = J GAA{£i,£2)fi£i)dei, (45) 

corresponding to the maximal possible eigenvalue Aq. 

As we apply a small magnetic field H, the percolation threshold C,c shifts away from 
its zero field value, However, according to our conjecture, the eigenvalue T remains 
constant. To first order in fisH and S^c, the correction to T is simply ST = (Ol^GlO), in 
complete analogy with first order perturbation theory in quantum mechanics. Thus the field 
induced correction to T is given by the following relation: 

1 + 7 J \7 dGAB[£2,£i) 7 oGaa[£i,£2) /\1 / 

(46) 

The condition 6T = may be rewritten in the form 

J f{e2m + j')6GAA{ei,e2)+-f{6GAB{ei,e2) + SGAB{e2,ei))]f{e,)de,de2 = 0. (47) 
This equation enables one to find the magnetic field induced correction to since 

GAAieue,) = G(^^,-'-^P-y (48) 

GAB{e^,e2) = G(^^,-'-^y (49) 

To first order in 6^c and (IbH, the condition 6T = (^) may be rewritten as 

£^12(0) \ f 2 H d{ei2 + e2i 



Jne2)f{e.)G'IC-'-^]\6Ul + ir-4 



dH 

From the definition (pO|) of EuiH), one finds 



deide2 = 0. (50) 



d 

— [euiH) + e2iiH)] = [1 + sign(£2 - ^i)] + /is [1 + sign(£i - €2)] = 2/iB. (51) 
Using Eqs. (0) and (§T|), we recover the result of the perturbative approach: 

= (1 + 7)^ T • (^2) 

To consider the case of arbitrarily strong magnetic fields, we change to the dimension- 
less variables Ai and A2 and, in doing so, reduce the integral eigenvalue equation to a 
dimensionless form which can then be solved numerically: 



[ hB{A2 



2) 
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f({l - Aum'Oil - (1 - A,,{x)fQ{^ - Ai2(x))\ {hA{A,)\ 

J lv7-'/'(l-A2i(x))'^(l- A2i(x)) 7~Ml-A2i(0))'^^^(l-A2i(0)) J [hBiA^) ) 

(53) 

Here, the eigenvalue is related to T by Eq. (pOl) . For the two-color problem, = 
Arf(x, 7), and the invariance conjecture is equivalent to the assumption that the product 
^d{x^ 7) l^ciXy 7)]"^^^ is independent of both the magnetic field and the ratio 7 of the two den- 
sities of states. Therefore using Eq. (PU]), the universal scaling function F{x) = ^c(T, H)/^^ 
can be expressed in terms of the eigenvalue A^, 



Ad(0,0) 
Ad(x,7) 



i/{d+i) 

(54) 



Formula (0) expresses F in terms of x = fJ'sH/T^c, instead of x = hbH/T^^. Noting that 
xlx = iciic = we find: 



x = x 



A.(0.0)1'"^«' 

,A^J ■ '''' 



Eqs. ([5^ ) and ( |55D determine F{x) in parametric form. 

In order to assess the validity of the invariant approach, we solved Eq. (p3|) numerically 
at the two values of 7 used in Section [V^, 7 = 1 and 1/2, over the entire relevant range of 
dimensionless magnetic fields. The resulting plots of F{x) are shown in Fig 2ab together 
with the two plots obtained via simulations. Let us first examine the two curves obtained 
for 7 = 1/2 (Fig 2a). It is clear that at weak fields the invariant method yields an accurate 
approximation of the universal scaling function F{x) at 7 = 1/2. In addition, we see that 
at 7 = 1/2 this approach gives the proper limiting value of F{x). At intermediate fields, the 
function F{x) obtained from the invariant method continues smoothly between its limiting 
forms, deviating from the corresponding values of F{x) obtained from simulations by at 
most 2%. 

Upon examination of the two curves in Fig 2b, it is equally clear that at 7 = 1 the 
invariant method accurately reproduces the scaling function F{x) at small fields. This 
method gives the proper limiting value of F{x) at 7 = 1 which, however, is higher than the 
corresponding simulation value by approximately 1.5%. This disagreement can be attributed 
to the systematic suppression of F{x) at large x due to certain finite size effects in the 
simulations run at 7 = 1 (see SectionjV^. The function F{x) obtained from the invariant 
method smoothly continues between its limiting forms, differing from the corresponding 
values of F{x) obtained from simulations by at most 6%. We believe that the magnitude of 
this discrepancy between the scaling functions, determined by simulations and the invariant 
method, at intermediate fields is considerably enhanced by the above mentioned finite size 
effects in simulations run at 7 = 1. 

We should point out that, for the values of F{x) obtained from the simulations, the 
standard error in F{x) is approximately 0.005. Thus at 7 = 1/2 there is a wide range 
of intermediate fields within which the difference between the numerically and analytically 
determined scaling functions is larger than the standard error. For 7 = 1 the two curves differ 
by more than the standard error for x > 0.05. Thus the invariant method does not suffice 
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to accurately determine the universal scaling function F{x). We see two possible reasons 
for this shortcoming: the first being that this procedure might not yield information about 
the infinite cluster and the second being the possibility that the average number of bonds 
per site, for sites belonging to the infinite cluster, is not an invariant of the dimensionless 
percolation problem. These questions are the subject of future work. 

VIII. CONCLUSION 

In this paper, we considered the problem of hopping magnetoresistance induced by Zee- 
man splitting. The problem has been reduced to the calculation of a universal function F{x) 
of a single dimensionless parameter x (which depends upon magnetic field and temperature). 
To find F{x), one has to solve certain dimensionless "two-color" percolation problem. We 
found the limiting behavior of F{x) analytically and obtained its full shape numerically for 
two values of the ratio Qa/qb = 1 and QaIqe = 1/2. In addition, we developed an approxi- 
mate method which enables us to solve the "two color" percolation problem analytically at 
any x. This approach gives a satisfactory extrapolation of the function F{x) between its 
two limiting forms. 

This theory can be applied to any system of localized electrons for which the width of 
the distribution function of the energies of the localized states is larger than the on-site 
Coulomb repulsion energy. Bulk amorphous silicon satisfies this condition . In the bulk 

= 3), there exists a range of temperatures (|^) within which Zeeman splitting makes the 
dominant contribution to the magnetoresistance. As was proposed in 0, this theory can 
then be used to probe the relative concentration of singly and doubly occupied sites which 
contribute to transport. 

This work was initiated by extensive communications with M. R. Beasley and D. Ephron. 
We are also indebted to A. S. loselevich and B. I. Shklovskii for helpful discussions. This work 
was supported by NSF Grant DMR-9423244. One of us (PC) acknowledges the Department 
of Education for fellowship support. 
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FIGURES 



FIG. 1. The universal function, F{x) for d = 2. The procedure for extracting F{x) and x 
from the simulation data gives the standard errors 0.005 for F{x) and less than 0.003 for x. Insets: 
normalized critical site concentration as a function of dimensionless magnetic field x- The standard 
error of ndx) / ^ci^) is 0.01 on average, (a): 7 = 1; for x > Xth, there is a systematic suppression 
of F{x) below its true value, F{xth) = 1-26; Inset: a linear fit to the data for x < 0.20 gave a 
slope of 1.47 ± 0.03 in good agreement with the theoretically predicted value, 3/2. (b):7 = 1/2; 
for X > Xth, the function F(x) attains its exact value, -F(xth) = 1-14; Inset: a linear fit to the data 
taken at X < 0.2 yielded a slope of 1.33 ±0.04 in good agreement with the value 4/3 following from 
the perturbation approach. 

FIG. 2. The approximate scaling function F(x) obtained via the invariant approach, for d = 2, 
shown with the corresponding curves extracted from simulations, (a): 7 = 1; the analytically 
derived scaling function reproduces the simulation results only at very small values of x. The 
difference in the limiting values of the two curves can be attributed to the finite size effect, described 
earlier, which is present in simulations run at 7 = 1. In addition , this finite size effect broadens the 
range of the scaling variable x, within which the two curves differ substantially, (b): 7 = 1/2; for 
X <^ Xth and x > Xth, the function F(x) obtained via the invariant approach is in good agreement 
with the numerically determined function. Within a large range of intermediate values of the 
scaling variable x, however, the two functions differ by more than 0.005, the standard error of the 
numerically determined F{x). 



16 



